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DYNAMICS  OF  IMPLODING  NEON  GAS  PUFF  PLASMAS  - 1 


I.  Introduction 


Recent  advances  in  the  technology  for  producing  gas  puff 
plasmas  has  made  them  very  attractive  as  high  brightness 
laboratory  x-ray  radiation  sources.  In  addition,  it  has 
become  evident  that  the  gas  puff  technology  can  be  employed 
to  investigate  a  number  of  interesting  x-ray  laser  schemes. 
In  fact,  Sandia  National  Laboratory  has  carried  out  a 
program  involving  the  implosion  of  a  gas  puff  plasma  onto  a 
low  density  foam  to  determine  the  feasibility  of  creating  a 
homogeneous,  uniform  gain  medium  along  the  central  core 
axis.1  The  preliminary  experimental  results  indicate  high 
radiative  conversion  efficiencies  from  relatively  clean, 
reproducible  implosions. 

In  this  investigation,  which  is  the  first  in  a  series  of 
reports,  we  seek  to  determine  a  better  understanding  of 
imploding  gas  puff  plasmas  for  a  varity  of  conditions  and 
configurations,  particularly  the  radiation  hydrodynamics 
with  and  without  a  central  core  plasma.  We  will  also 
explore  the  consequences  of  invoking  a  number  of 
approximations  and  their  impact  on  the  results.  For 
instance,  are  radiative  losses  of  sufficient  magnitude  to 
warrant  a  self-consistent  radiation  hydrodynamics  treatment? 
How  does  opacity  affect  the  overall  implosion  dynamics? 
Does  the  transport  of  energy  in  the  plasma  modify  the 
implosion  hydrodynamics?  What  detail  is  required  to 
accurately  model  the  radiation?  How  does  LTC  or  corona 
equilibrium  compare  with  CRE?  We  hope  to  shed  some  light  on 
these  and  other  issues  relevant  to  an  accurate  description 
and  understanding  of  imploding  gas  puff  plasmas. 

Manuacript  approved  Mar  8, 1988. 


II.  Theoretical  Model 


Since  the  implosion  of  a  cylindrical  gas  puff  plasma  can 
result  in  a  substantial  fraction  of  the  total  plasma  energy 
being  radiated  away,  the  pure  hydrodynamic  evolution  of  the 
plasma  may  be  modified.  Thus,  because  of  the  nonlinear 
nonlocal  interactions,  the  hydrodynamic  development,  atomic 
and  radiation  physics  of  the  plasma,  as  well  as  the 
transport  of  radiation,  must  be  calculated  self- 
consistently . 

Discussion  of  the  theoretical  model  can  be  separated  for 
convenience  as  follows:  (a)  hydrodynamics  and  thermal 
conduction  (b)  ionization  and  atomic  physics  and  (c) 
radiation  emission  and  transport. 


(A)  Hydrodynamics  and  Thermal  Conduction 

The  basic  hydrodynamic  variables  of  mass,  momentum,  and 

total  energy  are  transported  in  one  dimension  using  a 

numerical  scheme  with  a  sliding- zone  version  of  flux- 

2 

corrected  transport.  A  special  gridding  algorithm  is  used 
which  moves  zones  in  a  Lagrangian  fashion  and  adjusts  the 
mesh  in  order  to  resolve  steep  gradients  in  the  flow.  The 
hydrodynamic  equations  solved  are 
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where  p  is  mass  density,  u  is  velocity,  P  is  pressure,  ex  is 
total  energy  density,  era<j  is  the  rate  of  energy  loss  or 
gain  due  to  radiation,  h  is  the  thermal  conductivity,  and  N 
is  the  ion  density.  The  thermal  conduction  is  calculated 
implicitly,  using  an  iterative  Crank -Nicholson  scheme. 

Since  the  density  generally  did  not  exceed  solid  density 
in  this  study,  a  simple  equation  of  state  was  assumed,  viz. 

p  ■  I  (£t  "  i  pu2  -  CI)  '  (4) 

where  ej  is  the  potential  energy  due  to  ionization  and 
excitation.  (A  non- ideal  equation  of  state  taking  account  of 
ionization  energy  and  degeneracy  pressure  can  be  employed  in 
cases  where  the  density  exceeds  solid  density.)  A  single 
temperature  model  was  employed, 

kT  -  - - - 3-  ,  (5) 

(P/ mx)  ( 1+Z ) 

where  mx  is  ion  mass,  and  T  is  temperature.  The  ionization 
energy,  ej,  and  effective  charge,  Z  are  calculated  from  the 
ionization- radiation  equations  which  are  explained  below.  A 
single  temperature  assumption  is  valid  in  the  core  plasma, 
where  the  equilibration  time  is  of  the  order  of  picoseconds, 
and  it  is  adequate  in  the  stagnation  region,  where  the 
equilibration  time  can  be  of  the  order  of  nanoseconds.  In 
the  blowoff  plasma,  it  is  a  marginal  approximation,  but  the 
consequences  are  minor,  since  little  radiation  is  emitted 
from  this  region,  and  most  of  the  thermal  energy  is  carried 
by  the  electrons  in  the  blowoff. 

The  local  rate  of  change  of  energy  due  to  radiation 
transport,  era<j,  will  be  discussed  below. 
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(B)  Ionization  and  Atomic  Physics 

The  ionic  populations  in  the  plasma  are  determined  by  a 
set  of  atomic  rate  equations  of  the  form 

ar  '  J  Mii  f5  -  f  mij  £i  161 

where  fi  is  the  fractional  population  of  atomic  level  i,  and 
Wji  is  the  net  reaction  rate  describing  the  transition  from 
initial  state  j  to  final  state  i.  An  equation  of  this  type 
is  constructed  for  each  of  the  atomic  levels  included  in  the 
model . 

For  sufficiently  dense  plasmas,  the  effective  populating 
and  depopulating  rates  are  generally  fast  compared  with  the 
hydrodynamic  response.  Under  these  circumstances,  an 

equilibrium  assumption  can  be  justified,  which  involves 
dropping  the  explicit  time  dependence  in  Eq.  (6).  The 
plasma  is  then  said  to  be  in  collisional - radiative 
equilibrium  (CRE),  whereby  the  plasma  ionization  state 
responds  instantaneously  to  changes  in  hydrodynamic 

quantities. 

The  rate  coefficients  that  are  used  to  calculate  the 

populating  and  depopulating  rates,  Wji,  are  calculated  using 

various  atomic  calculational  methods.  The  processes 

included  in  this  calculation  and  the  methods  used  in 

calculating  the  corresponding  rate  coefficients  are 

4-10 

summarized  elsewhere. 

Once  the  set  of  rate  equations  (including  the  optical 
pumping  from  the  radiation  field)  has  been  solved  for  the 
level  populations  fi,  the  electron  density  can  be 

calculated, 
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where  Zi  is  the  ionic  charge  of  level  i  and  Ni  is  the  total 
ion  density. 


The  ionization  and  excitation  energy  can  also  be 
calculated  by 
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wh<jre  Xi  is  the  energy  of  level  i,  measured  from  the  ground 
state  of  the  neutral  atom. 


For  the  simulations  presented  below,  the  atomic  model 
for  neon  contains  27  atomic  levels  and  13  emission  lines. 


(C)  Radiation  Emission  and  Transport 

Radiation  emission  from  and  absorption  by  a  plasma  are 
dependent  on  the  local  atomic  level  population  densities. 
Except  for  optically  thin  plasmas,  however,  the  level 
populations  depend  on  the  radiation  field,  since  optical 
pumping  via  photoionization  and  photoexcitation  can  produce 
significant  population  redistribution.  Thus,  the  ionization 
and  radiation  transport  processes  are  strongly  coupled  and 
must  be  solved  self -consistently.  In  this  model,  an 
iterative  procedure^  is  used,  where  level  populations  are 
calculated  using  the  radiation  field  from  the  previous 
iteration,  then  using  these  populations  to  calculate  a  new 
radiation  field  and  recalculating  populations  until 
convergence  is  reached. 
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A  probabilistic  radiation  transport  scheme  ' 
employed,  which  forms  local  angle  and  frequency  averaged 
escape  probabilities  for  each  emission  line  and  for  each 
bound-free  process.  Free-free  radiation  is  treated  with  a 
multi  frequency  transport  formalism.  The  radiation  transport 
and  emission  spectra  are  calculated  from  these  escape 
probabilities.  The  method  can  treat  comprehensive  atomic 
models  and  provides  good  overall  energetics,  but  cannot 
calculate  accurately  certain  spectral  details  and  lines  with 
very  high  optical  depths. 

Inner-shell  opacities  are  included  in  the  model,  since 
these  processes  are  very  important  in  the  cool,  dense  plasma 
regions.  Inner-shell  photoionization  cross  sections  for  the 
neutral  element  are  taken  from  the  fits  by  Biggs  and 
Lighthill,14  and  the  positions  of  the  ionization-dependent 
absorption  edges  are  taken  from  the  Hartree-Fock 
calculations  of  dementi  and  Roetti.*^ 

The  local  rate  of  energy  change  in  zone  j,  due  to 
radiation  transport  is  given  by 
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where  Fpfc  is  the  rate  of  energy  loss  in  zone  k  due  to  a 
discrete  radiative  process  (or  frequency  group)  P,  and  Cpkj 
is  the  radiative  coupling  of  zone  k  to  zone  j  for  that 
process.  The  couplings  are  functions  of  opacity,  integated 
over  process  and  photon  path.  In  the  probabilistic  model,  a 
matrix  of  couplings  must  be  computed  for  each  bound-bound, 
bound-free  and  free-free  process.  In  this  way,  the  net 
cooling  and  heating  by  radiation  emission  and  absorption 
between  the  various  zones  of  the  plasma  is  accurately  taken 
into  account. 
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III.  Results 

(A)  Puff-Gas  Implosion  with  a  Central  Core  Plasma 

Simulations  were  performed  for  a  cylindrical  annular 
neon  puff -gas  of  density  5X10’^  g/cm^,  with  inner  radiu-s 
0.55  cm  and  outer  radius  1.95  cm,  imploding  radially  at  a 
velocity  of  3X10^  cm/sec.  Results  were  obtained  with  and 
without  a  central  neon  core  plasma  of  density  5X10"^  g/cm^ 
and  radius  0.10  cm.  The  initial  configuration  with  the  core 
plasma  is  shown  in  Fig.  1.  The  temperature  of  the  puff -gas 
was  taken  to  be  about  5  eV  initially,  and  the  core  plasma 
temperature  was  about  0.04  eV.  A  tenuous  background  plasma 
was  placed  between  the  puff -gas  and  core  plasmas  of  density 
5X10"^  g/cm*.  These  initial  conditions  were  chosen  to 
correspond  approximately  to  experiments  being  conducted  at 
Sandia  National  Laboratory.* 

Results  of  simulations  with  a  central  core  plasma  will 
be  discussed  first.  For  the  first  15  nanoseconds,  the  puff- 
gas  plasma  essentially  coasts  radially  inward.  The  forward 
edge  of  the  puff  gas  is  heated  to  a  few  tens  of  eV  through 
accretion  of  background  plasma.  At  about  15  nanoseconds, 
the  puff -gas  makes  full  contact  with  the  core,  and 
conversion  of  kinetic  energy  to  thermal  energy  at  the 
interface  creates  a  large  overpressure.  The  plasma 
temperature  quickly  exceeds  100  eV  in  the  contact  region, 
and  a  narrow  region  of  intense  net  radiative  emission 
centered  at  the  interface  is  formed. 

At  40  nanoseconds  (Fig.  2),  the  peak  temperature  near 
the  interface  exceeds  300  eV,  and  the  overpressure  has 
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caused  shocks  to  propagate  radially  inward  in  the  core 
plasma  and  outward  in  the  puff  gas.  Heating  from  the  soft 
x-rays  produced  in  the  interface  region  is  evident  in  the 
dense  core,  and  radiation  heating  is  also  taking  place  deep 
in  the  puff  gas.  The  effects  of  radiation  transport  will  be 
more  clearly  seen  when  comparison  is  made  with  a  radiation¬ 
less  simulation.  Thermal  conduction  creates  a  nearly 
isothermal  region  near  the  leading  edge  of  the  puff  gas,  but 
it  is  ineffective  in  transporting  energy  into  the  core. 

By  50  nanoseconds  (Fig.  3),  the  inward  propagating  shock 
has  reached  the  origin.  Although  the  plasma  approaches 
solid  density  on  axis,  the  temperature  there  is  still  only  a 
few  tens  of  eV.  This  region  of  very  dense  plasma  is 
radiating  strongly  in  the  L-shell,  but  its  volume  is  small. 
A  broader  annulus  of  hot  plasma,  situated  outside  the  core 
plasma  is  producing  most  of  the  K-  and  L-shell  radiation. 
The  emission  spectrum  at  50  nanoseconds  is  shown  in  Fig.  4. 
The  spectral  power  is  about  evenly  divided  between  continuum 
and  line  radiation  at  this  time;  almost  half  of  the  total 
power  is  carried  by  the  K-shell  emission  lines  Ne  x  2p-ls 
and  Ne  IX  2p-ls. 

At  70  nanoseconds  (Fig.  5)  the  entire  plasma  is  moving 
radially  outwards.  The  core  plasma  remains  cool,  since  the 
density  is  still  high  enough  to  retard  thermal  conduction. 
Furthermore,  the  radiative  emission  of  this  region  is  in 
approximate  balance  with  the  radiative  absorption. 

The  energy  history  of  the  plasma  throughout  the  first 
120  nanoseconds  is  shown  in  Fig.  6.  Initially,  about  98%  of 
the  total  energy  is  kinetic;  the  remainder  is  in  thermal  and 
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ionization  energies.  The  energy  lost  to  radiation  begins  to 
increase  sharply  just  before  20  nanoseconds,  as  the  puff  gas 
comes  into  contact  with  the  central  core.  Thermal  and 
ionization  energies  also  begin  to  increase  sharply  at  this 
time.  Kinetic  energy  goes  through  a  well  defined  minimum  at 
about  60  nanoseconds,  at  which  time  almost  half  of  the  total 
plasma  energy  has  been  radiated  away.  Although  peak 
compression  on  axis  occurs  at  about  50  nanoseconds, 
substantial  portions  of  the  puff  gas  plasma  continue  to  move 
radially  inward  until  about  60  nanoseconds.  By  120 
nanoseconds  about  two-thirds  of  the  total  energy  has  been 
converted  to  radiation. 

A  similar  calculation  was  performed  without  radiation 
effects,  but  with  detailed  atomic  physics,  so  that  the 
ionization  and  excitation  energies  and  electron  number 
densities  would  be  self -consistently  calculated.  Although 
radiative  cooling  and  optical  pumping  were  turned  off, 
optically  thin  radiation  was  calculated  for  comparison  with 
the  simulation  discussed  above. 

At  40  nanoseconds  (Fig.  7),  the  profiles  are 
qualitatively  very  similar  to  those  with  radiation  effects. 
However,  the  role  of  radiation  in  transporting  energy  is 
clearly  evident.  The  dense  core  is  still  very  cold  compared 
with  the  complete  simulation  with  radiation  transport. 
Also,  the  temperatures  in  the  puff-gas  plasma  differ  by  a 
factor  of  about  three.  The  heating  of  the  puff-gas  out  to 
about  0.8  cm  evident  in  Fig.  2  is  the  direct  result  of 
radiation  transport.  Temperature  and  density  gradients  in 
the  inward  propagating  shock  are  steeper  in  the 
radiationless  case;  however,  assembly  occurs  at  about  the 
same  time. 
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HYDRO  PROFILES 
20.00  nsec 


TEMPERATURE 
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DISPIACEMENT 

Profiles  at  20  nanoseconds  for  radiation-hydrodynamics  calculation 
without  a  central  core.  Temperature  and  density  at  the  origin  are 
near  peak  values;  as  plasma  continues  to  flow  inward,  the  volume 
of  hot  moderately  dense  plasma  will  increase. 


Fig.  9. 
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thin  rod 
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Fig.  8.  Energy  partition  for  puff  gas  implosion  with  a  central  core,  but 
neglecting  the  effects  of  radiation.  Now  the  total  energy  is 
synonymous  with  the  plasma  energy.  Optically  thin  radiation  is 
plotted  for  purposes  of  comparison,  but  it  does  not  remove  energy 
from  the  plasma.  The  relative  minimun  in  kinetic  energy  occurs 
somewhat  earlier,  and  the  energies  are  substantially  different 
from  those  in  Fig.  6  after  30  nanoseconds. 
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Profiles  at  HO  nanoseconds  with  a  central  core  neglecting  the 
effects  of  radiation.  The  core  plasma  remains  very  cold  and  there 
is  no  evidence  of  heating  deep  in  the  puff-gas. 
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Energy  partition  for  puff  gas  implosion  with  a  central  core. 
Total  energy,  which  remains  constant  in  time,  is  the  sum  of  the 
energy  lost  to  radiation  and  the  energy  which  remains  in  the 
plasma.  The  plasma  energy  is  the  sum  of  the  kinetic,  ionization 
and  thermal  energies.  The  energies  (log10  energy  in  ergs)  are 
given  as  a  function  of  time  (in  nanoseconds). 


displacement 


Profiles  at  70  nanoseconds  for  radiation-hydrodynamic  calculation 
with  a  central  core.  Expansion  phase:  the  entire  plasma  is  now 
moving  radially  outwards.  The  core  plasma  remains  cool. 


.01  0.1  1.0  10.0 


.  k.  Emission  spectrum  at  50  nanoseconds.  Spectral  intensity  (from 

10' 5  to  102  ergs/3ec-cm-keV)  is  plotted  as  a  function  of  photon 
energy  (from  10  eV  to  10  keV).  K-shell  radiation  dominates  at 
this  time  with  most  of  the  power  in  the  indicated  H  and  He-like 
emission  lines. 


hydro  profiles 
50.00  NSEC. 


displacement 


Fig.  3-  Profiles  at  50  nanoseconds  for  radiation-hydrodynamics  calculation 
with  a  central  core.  Assembly  phase:  inward  propagating  shock 
reaches  axis;  peak  density  approaches  solid  density.  Bulk  of 
plasma  is  still  moving  radially  inwards. 
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Profiles  at  40  nanoseconds  for  radiation-hydrodynamics  calculation 
with  a  central  core.  Inward  and  outward  propagating  shocks  have 
been  formed.  Radiation  transport  has  heated  core  plasma  and  puff 
gas  out  to  about  0.8  cm.  Radial  displacement  is  given  in  cm., 
log10  density  in  g/cnH,  iog10  temperature  in  eV.log^  pressure  in 
ergs/cnH  x  10“  2.  log10  velocity  in  cm/nsec  x  10-2  and  log1Q 
e  .  in  ergs/cm^-nsec  x  10”  . 


Fig.  1.  Initial  Configuration:  Puff  gas  plasma  is  imploding  at  0.03 
cm/nsec  onto  a  low  density  core. 
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hydrodynamic  response,  a  self-consistent  calculation, 
linking  radiation  and  fluid  transport  with  the  ionization 
dynamics,  must  be  performed.  The  "hydrodynamic" 
calculation,  which  was  compared  to  the  self-consistent 
simulation,  neglected  radiation  but  included  ionization 
dynamics.  At  certain  times  and  regions,  the  local  energy  of 
ionization  and  excitation  substantially  exceeded  the  local 
thermal  energies.  A  pure  "hydrodynamic"  calculation  with 
simplified  prescriptions  for  ionization  energy  and  average 
ionization  state  would  produce  a  substantially  different 
plasma  evolution,  leading  to  erroneous  conclusions. 

Since  radiation  plays  such  a  crucial  role,  it  must  be 
generated  and  transported  accurately.  As  we  have  shown, 
optically  thin  radiation  ignores  the  importance  of  transport 
of  radiative  energy  and  grossly  overestimates  losses.  LTE 
or  corona  equilibrium  models  for  the  level  populations  will 
also  lead  to  inaccurate  energetics.  We  found  that  the 
atomic  populations  in  the  regions  of  peak  radiative  emission 
were  substantially  different  from  corresponding  LTE  or 
corona  equilibrium  populations,  in  some  cases  by  orders  of 
magnitude . 

Accurate  radiation  energetics  demands  a  sufficiently 
detailed  atomic  model.  The  simple  neon  model  employed  here 
underestimates  the  radiation  losses  by  about  a  factor  of 
two.  Thus,  we  plan  to  perform  similar  calculations  with  a 
full  neon  model  to  assess  and  determine  its  influence  on  the 
plasma  evolution,  radiative  yields,  and  spectrum. 


radiation  losses  under  these  circumstances.  Ultimately,  the 
optically  thin  treatment  overestimates  energy  loss  by  almost 
a  factor  of  two. 

The  simulations  described  above  were  carried  out  with  a 
relatively  simple  atomic  model  for  neon.  How  much  accuracy 
has  been  given  up  by  using  a  reduced  model?  The 
hydrodynamic  profiles  of  Fig.  9  were  post-processed  using  a 
full  neon  K-  and  L- shell  model,  and  the  resulting  spectrum 
is  shown  in  Fig.  14.  Not  only  are  there  more  emission 
lines,  but  the  radiative  power  in  some  of  the  most  important 
lines  is  substantially  increased,  as  is  the  continuum 
radiation.  The  ionization  state  of  the  plasma  is  actually 
increased,  reflecting  the  increased  probability  of 
ionization  from  the  excited  levels  (which  are  included). 
The  total  radiated  power  is  slightly  more  than  a  factor  of 
two  larger  than  that  calculated  with  the  simple  model. 

IV.  Conclusions 


Puff -gas  implosions  show  promise  as  a  means  of 
efficiently  converting  kinetic  energy  to  radiation.  Higher 
plasma  densities  and  total  radiative  yields  can  be  achieved 
through  the  use  of  a  central  core  plasma.  In  the  case  we 
studied  using  a  core,  more  than  two-thirds  of  the  total 
plasma  energy  was  converted  to  K-  and  L-shell  radiation. 
Without  a  core  plasma,  total  radiative  yield  is  smaller,  but 
the  fraction  of  K-shell  radiation  is  larger,  due  to  the 
reduced  plasma  mass. 

Because  the  radiative  energy  loss  and  the  transport  of 
radiation  in  the  plasma  substantially  modify  the 
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emission  spectrum  at  20  nanoseconds  is  given  in  Fig.  10. 
More  than  75%  of  the  spectral  power  is  carried  by  the  R- 
shell  lines.  Plasma  continues  to  flow  inward  and  stagnate, 
and,  although  the  peak  density  falls  somewhat,  by  40 
nanoseconds,  a  substantially  larger  volume  of  plasma  is 
radiating.  Since  the  temperature  remains  about  the  same, 
the  emission  spectrum  at  40  nanoseconds  is  qualitatively 
very  similar,  but  shifted  upward  in  intensity  by  about  an 
order  of  magnitude  {Fig.  11). 

The  energy  history  of  the  implosion  is  given  in  Fig.  12. 
It  is  similar  to  that  with  a  central  core  plasma  (Fig.  6), 
especially  during  the  first  50  nanoseconds.  In  the  present 
case,  the  radiated  energy  rises  even  more  rapidly,  and 
actually  exceeds  the  yield  with  a  central  core  until  about 
50  nanoseconds.  Ultimately,  slightly  more  than  half  of  the 
total  plasma  energy  is  radiated  away.  The  plasma  kinetic 
energy  exhibits  a  more  pronounced  minimum  at  a  slightly 
later  time.  Thus,  the  fraction  of  the  total  energy  radiated 
away  is  significantly  larger  with  a  central  plasma  core,  but 
the  additional  mass  results  in  a  smaller  average 
temperature,  and  the  fraction  of  K-shell  radiation  is 
substantially  smaller. 

A  simulation  was  performed  for  the  case  without  a  core 
neglecting  radiation  effects.  The  energy  history  for  this 
calculation  is  shown  in  Fig.  13.  Because  of  the  absence  of 
radiative  cooling,  thermal  energies  are  higher,  and  the 
kinetic  energy  minimum  occurs  sooner.  Optically  thin 
radiation,  which  was  calculated  but  not  used  to  cool  the 
plasma,  is  observed  to  rise  rapidly,  even  before  ten 
nanoseconds.  Opacity  effects  would  substantially  reduce 
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A  corresponding  energy  history  of  this  simulation  is 
shown  in  Fig.  8.  The  early  history  is  similar  to  that  shown 
in  Fig.  6.  By  30  nanoseconds,  however,  substantially  more 
energy  resides  in  thermal  and  ionization  energies.  The 
relative  minimum  in  kinetic  energy  occurs  earlier  and  is  not 
as  sharp.  The  difference  in  timing  can  be  attributed  to  the 
role  of  radiation  in  reducing  the  overpressure  in  the 
interface  region.  As  can  be  seen  from  the  optically  thin 
radiation  curve,  radiative  cooling  can  be  greatly 
overestimated  in  problems  of  this  type  if  opacity  effects 
are  not  included. 


Puff-Gas  Implosion  Without  a  Core  Plasma 


A  series  of  simulations  without  a  core  plasma  were 
performed,  assuming  the  same  puff  gas  density  of  5X10’^ 
g/cm^  and  temperature  of  5  eV.  As  before,  the  inner  and 
outer  radii  of  the  puff  gas  were  taken  to  be  0.55  cm  and 
1.95  cm,  respectively,  with  a  radial  velocity  of  3X107 
cm/sec.  A  background  plasma  of  density  5X10' 7  g/cm^  filled 
the  central  void. 

During  the  coasting  phase,  lasting  about  18  nanoseconds, 
the  puff  gas  moves  radially  inward,  accreting  background 
plasma  and  warming  to  a  few  tens  of  eV  at  its  forward  edge. 
Rapid  heating  and  compression  take  place  as  the  puff  gas 
reaches  the  origin.  Figure  9  shows  the  situation  at  20 
nanoseconds.  Peak  density  is  about  1.3  X  10"^  g/cm^  and 
temperature  near  the  axis  is  slightly  above  300  eV.  The 
hot,  moderately  dense  plasma  near  the  origin  is  radiating 
strongly  in  the  K-shell,  but  its  volume  is  small.  The 
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Fig.  12.  Energy  partition  for  puff  gas  implosion  without  a  central  core. 

Although  radiation  energy  increases  more  rapidly  at  early  times 
without  a  core,  only  about  half  of  the  total  energy  is  radiated 
away  by  100  nanoseconds,  compared  with  about  two-thirds  in  the 
case  with  a  core. 
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Fig.  14.  Emission  spectrum  at  40  nanoseconds  with  complete  atomic  model 

for  neon.  Comparison  should  be  made  with  Fig.  11.  Substantially 
increased  radiative  power  in  both  emission  lines  and  continuum. 
Ionization  of  plasma  is  increased  due  to  inclusion  of  more  excited 
levels. 
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